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ABSTRACT 


The Apollo Lunar Program utilized efficient, i.e., Earth-return, transearth 
trajectories which employed parking orbits in order to minimize energy requirements. 
This thesis concentrates on a different type of transearlh trajectory. These are direct- 
ascent, hyperbolic trajectories which omit the parking orbits in order to achieve short 
flight times to and from a future lunar base. The object of the thesis is the development 
of a three-dimensional transearlh trajectory model and associated computer program for 
exploring trade-offs between flight-lime and energy, given various mission constraints. 
The program also targets the Moon with a hyperbolic trajectory, which can with a time- 
reversed trajectory; be used for targeting Earth impact points. The first-order model is 
based on an Earth-centered conic and a massless spherical Moon, using MathCAD 
version 3.0. This model is intended as the basis for future patched-conic formulation for 
the design of fast Earth-return trajectories. Applications include placing nuclear- 
deterrent arsenals on the Moon, various space support related activities and finally 
protection against Earth-threatening asteroids and comets using lunar bases. 


C.,: 

TABLE OF CONTENTS 


I. INTRODUCTION . 1 

A. BACKGROUND . 1 

B. OBJECTIVES AND LIMITATIONS. 2 

II. ANALYSIS APPROACH. 3 

A. TRANSLUNAR TRAJECTORY DESIGN . 4 

1. Translunar Trajectory Types . 5 

2. Non-Inertial Coordinate Frames . 5 

B. SOLUTION FOR POINT-MOON INTERCEPT. 7 

1. Three Dimensional Point-Moon Model . 8 

2. Three Dimensional Considerations. 8 

C. SOLUTION FOR MASSLESS SPHERICAL MOON INTERCEPT . . 12 

1. Iterative Technique for Lunar Surface Impact . 12 

2. Selenographic Coordinates . 14 

a. Euler angles. 14 

b. Lunar Impact. 14 

D. PATCHED CONIC APPROXIMATION. 17 

III. APPLICATIONS. 19 


IV 


















A. STRATEGIC PEACE INITIATIVE CONCEPT . 19 

1. Background Information. 19 

2. SPI Architecture. 21 

B. SPI TECHNICAL FEASIBILITY . 22 

C. SPACE SUPPORT ACTIVITIES . 24 

D. SUMMARY. 24 

IV. CONCLUSIONS AND RECOMMENDATIONS. 26 

A. FINDINGS. 26 

B. RECOMMENDATIONS. 26 

APPENDIX A. PATCHED-CONIC BASED ON INERTIAL FRAMES .27 

APPENDIX B. LUNAR TRAJECTORY PROGRAM . 28 

LIST OF REFERENCES.55 

INITIAL DISTRIBUTION LIST. 56 


v 














TABLE OF SYMBOLS AND ABBREVIATIONS 


This listing of symbols and abbreviations is to aid the reader in clarifying the 
multitude of symbols and abbreviations used in this thesis. They are in alphabetical order. 

a E orb : input value for posigrade or retrograde orbit 

a E : spacecraft right ascension at Earth launch 

a M : spacecraft right ascension at Moon intercept 

a M1 : represents cosine of a M 

or M2 : represents sine of a M 

a EM : represents the difference between a M and a E 
a EMi : represents the cosine of o; EM 

a EM 2 : represents the sine of a EM 

ag, : represents the difference between the translunar trajectory longitude-of- 

ascending-node and spacecraft right ascension at Earth launch 
: represents the cosine of 
a EU 2 : represents the sine of a m 

a, : represents an intermediate value to calculate the spacecraft right 

ascension on its trajectory 

a 0 : represents spacecraft right ascension on its trajectory 

a : represents the selenographic longitude for lunar surface impact 

/3 e : represents flight path azimuth at launch 
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Pe\ 

: represents 

Pe2 

: represents 

/^M 

: represents 

Pm 

: represents 

Pm2 

: represents 

Po 

: represents 

P\ 

: represents 

P2 

: represents 

C E 

: represents 

Cei 

: represents 

C E 2 

: represents 

C M 

: represents 

C 0 

: represents 

Ye 

: represents 

Yei 

: represents 

Ye2 

: represents 

Ym 

: represents 

Ymi 

: represents 

Y M2 

: represents 

Yo 

: represents 

Yi 

: represents 

Y 2 

: represents 

D 

: represents 


cosine of /3, 
sine of jB b 

ilight path azimuth at moon intercept 
cosine of /3 M 
sine of /8 M 

[light path azimuth of the spacecraft on its trajectory 
cosine of /3 0 
sine of /3 n 

the sweep angle at launch 
the cosine of C E 
the sine of C E 

the sweep angle at moon intercept 

the sweep angle of the spacecraft on its trajectory 

Ilight path angle at launch 

the cosine of Ye 

the sine of Ye 

Ilight path angle at moon intercept 
the cosine of y m 
the sine of y.m 

the (light path angle of the spacecraft on its trajectory 

the cosine of Yo 

the sine of y 0 

the deflection angle 


VII 


D.mo 

5 m 

5o 

5 

Aa 

Aa M 

A® Ml 

Aa 0 
Aoq 
A a 2 
A p0 

ArO 


eccen 


eccen„ 


energy min 


represents the directrix-to-focus or axial offset of the asymptotic from 
the focus ; it is a variable valid for eccentricities greater than one 
represents the distance between the spacecraft and the lunar center at 
time t 0 

represents the spacecraft declination at earth launch point 
represents spacecraft declination at moon intercept 
represents the spacecraft declination on its trajectory 
represents the selenographic latitude for surface impact on the moon 
see Figure 8 

represents intermediate value to calculate spacecraft right ascension on 
its trajectory 

represents the cosine of Aa M 
represents the sine of Aa M 

represents an intermediate value for the calculation of a 0 
represents the cosine of Aa 0 
represents the sine of Aa 0 

represents cosine of X s : see iterative process on page 12 
represents the difference between point-moon intercept flight time 
relative to perigee passage and the spacecraft flight time on its 
trajectory 

represents eccentricity as an input value 

represents the minimum eccentricity input value allowed 

represents the minimum specific energy required for launch 


E hI : rcprcscnls ihc hyperbolic eccentric anomaly for earth launch 

E h2 : represents the hyperbolic eccentric anomaly for point-moon intercept 

E hJ : represents the hyperbolic eccentric anomaly for the spacecraft on its 

trajectory 

E shl : represents the hyperbolic sine of E h , 

E ih2 : represents the hyperbolic sine of E h2 

E sh3 : represents the hyperbolic sine of E h , 

f M : represents the input value for the true anomaly of the spacecraft at 

moon intercept 

f E : represents the true anomaly of the spacecraft at earth launch 

^EM : represents the difference of the true anomaly of the spacecraft at moon 

intercept and at earth launch 

I'm max : represents a constraint value on f M if required to be used 

f M_MAxi : represents the cosine of f M MAX 

f M.MAX 2 : represents the sine of f M MAX 

^Ml : represents the cosine of f M 

^M2 : represents the sine of f M 

f 0 : represents the true anomaly of the spacecraft on its trajectory 

f\io : represents the cosine of f fl 

0 : represents the first Euler angle 

h E : represents the geocentric altitude of launch point 

i M : represents the inclination of the moon’s orbit plane 

i : represents the translunar trajectory inclination 


IX 


i, : represents the cosine of i 

i 2 : represents the sine of i 

k E : represents the earth’s gravitational constant 

k M : represents the moon’s gravitational constant 

r E : represents the earth’s mean radius 

r M : represents the moon’s mean radius 

R E1 : represents the spacecraft radial distance at earth launch 

R M1 : represents the spacecraft radial distance at moon intercept 

r : represents the radial distance ratio 

R a : represents apogee distance : valid only for eccentricities less than 1 

R p : represents the perigee distance 

R 0 : represents the input value for the position of the spacecraft on its 

trajectory 

R m : represents the radial velocity at point-moon intercept 

R x : represents an intermediate value to calculate f 0 

Pem : represents the Earth-Moon mean distance 

sigrifjjM : represents an intermediate value calculated to insure correct use of sign 
when calculating a EM 

SfEM : represents the input value of +1 or - 1 to insure correct sign usage in 
calculating a EM 

T p : represents a circular orbit period for radius R p 

t E : represents Earth launch flight time relative to perigee passage 

t E1 : represents an intermediate value used to calculate t E 
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t M : represents the point-moon intercept llight time relative to perigee 

passage 

l M1 : represents an intermediate value used to calculate t M 

t EM : represents the point-moon translunar trajectory llight time 

t M0 : represents tin intermediate value used to calculate t 0 

t 0 : represents the (light lime of the spacecraft on its trajectory 

V p : represents the perigee velocity 

V E : represents the spacecraft velocity at launch 

V M : represents the spacecraft velocity at moon intercept 

V 0 : represents the spacecraft velocity on its trajectory 

4> m : represents the input value of the lunar sweep angle at moon intercept 

c f> m : represents the cosine of <£ M 

< f> M2 '■ represents the sine of <£ M 

$ M0 : represents the difference between the lunar sweep angle at moon 

intercept and the moon's mean orbital rotation rate multiplied by A^ 
<fi : represents the first Euler angle 

X M : represents the X-coordinate for the rectangular coordinates of the moon 

at point-moon intercept 

X M - : represents the X-component for calculating radial velocity at point-moon 
intercept 

X 0 : represents the X-coordinate for the rectangular coordinates of the 

spacecraft on its trajectory' 

X M0 : represents the X-coordinate for the center of the moon 


XI 


X : represents the difference between X M and X M0 

X| un : represents the X-coordinate for the rectangular selenographic 
coordinates 

Y m : represents the Y-coordinate for the rectangular coordinates of the moon 

at point-moon intercept 

Y m . : represents the Y-component for calculating radial velocity at point-moon 
intercept 

Y 0 : represents the Y-coordinate for the rectangular coordinates of the 

spacecraft on its trajectory' 

Y mo : represents the Y-coordinate for the center of the moon 

Y : represents the difference between Y M and Y M0 

Y lun : represents the Y-coordinate for the rectangular selenographic 
coordinates 

’P M : represents the third Euler angle 

Z M : represents the Z-coordinate for the rectangular coordinates of the moon 

at point-moon intercept 

Z M : represents the Z-component for calculating radial velocity at point-moon 

intercept 

Z M0 : represents the Z-coordinate for the center of the moon 

Z : represents the difference between Z M and Z M0 

Z, un : represents the Z-coordinate for the rectangular selenographic 
coordinates 

Cl : represents the translunar trajectory longitude-of-the-ascending-node 


xii 


u>, : represents ihe Earth's angular rotational rate 

u) M : represents the Moon’s mean orbital rotation rate 

u> : translunar trajectory argument-of-perigee 
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I. INTRODUCTION 


This chapter provides a general description of this thesis and is divided into two 
sections. The first section presents background information pertaining to translunar 
trajectory design. The second section delineates the objectives and limitations of this 
project. 

A. BACKGROUND 

While the impending mission of establishing a permanent lunar base stands in our 
future, it is appropriate to study this type of mission by investigating various translunar 
trajectories. The trade-off between flight lime and energy required must be determined 
prior to considering the type of translunar trajectory to design. The energy required for 
the hyperbolic trajectory can then be compared to that of the Apollo mission to 
determine the increase in energy and thereby cost increases. 

The Apollo Lunar Program of the 1960s utilized energy-efficient translunar 
trajectories. This thesis will analyze a different class of trajectories named fast Earth- 
return trajectories. These trajectories utilize a direct-ascent from the lunar surface and 
omit the lunar parking orbit used in the Apollo missions. These direct-ascent trajectories 
differ considerably from translunar trajectories previously considered and discussed in this 
thesis. (Wadsworth, 1991) 
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B. OBJECTIVES AND LIMITATIONS 


The objective of this thesis is to develop a trajectory model and associated 
computer program for systems studies of fast Earth-return trajectories from a lunar base. 
The scope of the project involves a step process which evolves in complexity, as the 
design moves from a massless point-Moon, to a massless spherical Moon intercept 
solution. Finally, the groundwork is laid for a patched-conic model which culminates in 
a solution for intercept on a Moon with finite lunar mass. A novel feature is the use of 
a time-reversed trajectory for targeting Earth impact points. In other words, the lunar 
impact point serves as the lunar launch point for a time-reversed trajectory. The 
complete patched conic approximation, as discussed in Chapter II, is not developed in the 
program due to scope of work and time constraints. 
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II. ANALYSIS APPROACH 


This chapter sets the foundation of the thesis by delineating a How chart of tasks 
that are required to be accomplished to reach the objective of the thesis. A flow chart 
represented in Figure 1, summarizes the work to be accomplished. The last block of the 
diagram represents further research work. 



Figure 1 . Translunar Trajectory Design Flowchart. 

Precision translunar/transearth trajectory design requires time-consuming, step-by- 
step numerical solution of the exact differential equations of motion, including lunar and 
solar ephemerides. For preliminary design, the computation lime can be reduced a 
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hundred-fold or more by employing a closed-form "patched-conic" approximation to the 
precision trajectory. Since the patched-conic is a composite of the two conics, it can be 
formulated explicitly in terms of elementary functions of the key parameters. This 
provides better insight into design trade-offs than numerical computation. The patched 
conic approximates the gravitational effects of the Earth-Moon system. It ignores the 
small perturbation due to the Sun. 

Two approaches were considered for developing the patched-conic model. Both 
approaches employ Earth-centered and Moon-centered coordinate frames with respect 
to which the conic trajectories are formulated. Although the conic trajectories are planar 
with respect to the Earth-centered and Moon-centered frames, they form three- 
dimensional curves as viewed in an inertial barycentric frame. In the first approach, the 
coordinate frames are non-inertial, being attached to the Earth and Moon which travel 
in (idealized) circular orbits (with respect to the Earth-Moon inertial frame). The second 
approach employs non-rotating, inertial Earth-centered and Moon-centered coordinate 
frames. These frames are attached, respectively, to a fictitious Earth and fictitious Moon 
which travel at constant velocity with respect to the barycentric frame. The implications 
of the second approach are briefly stated in Appendix A. The advantage of the first 
approach, pursued in this thesis research, is described in section A. (Wadsworth, 1991) 

A. TRANSLUNAR TRAJECTORY DESIGN 

The first consideration given to designing a translunar trajectory is to decide on 
what translunar trajectory types will be used. Since this thesis concentrates on fast Earth- 
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return trajectories, eccentricities greater than one will be utilized. The reason for using 
non-inertial coordinate frames is discussed in section A.2. 

1. Translunar Trajectory Types 

Figure 2 illustrates all possible trajectory types that can be considered for 
Earth-Moon trajectory' design. The radius of the Earth (R,), radius of the Moon (R M ), 
spacecraft true anomaly at launch (f E ), and spacecraft true anomaly at moon intercept 
(f M ) are projected in Figure 2. 

Since only trajectories with eccentricities greater than 1.0 are being 
considered, only type I and III hyperbolic trajectories will be used in thesis research (see 
Table I). 

2. Non-inertial Coordinate Frames 

The Earth-centered and Moon-centered coordinate frames are non-inertial 
frames, since they are in a state of gravitational free-fall about the Earth-Moon 
barycenler (center of gravity). As a consequence, the composite palched-conic provides 
a better approximation to the restricted, three-body spacecraft motion than would be 
suggested by simply joining two conics. 

Within the Moon’s region-of-influence, the terrestrial perturbing acceleration 
of the spacecraft is nearly the same as the centripetal acceleration of the Moon-centered 
coordinate frame. In this frame, these accelerations nearly cancel; this can be deduced 
from the facts that the distances of the spacecraft to the Earth and the Moon to the 
barycenler are comparable and the mass of the Earth is comparable to the equivalent 
mass at the barycenler. For the Earth-centered phase, there is a similar benefit over 
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Figure 2. Equatorial plane projection. 


TRANSLUNAR TRAJECTORY TYPES: 


TYPE I posigrade : arc ab 
TYPE I retrograde : arc dc 
TYPE II posigrade . arc abc ( e < 1 ) 
TYPE II retrograde : arc deb ( e < 1 ) 
TYPE III posigrade : arc dab 
TYPE III retrograde : arc adc 
TYPE IV posigrade : arc dabc ( e < 1 ) 
TYPE IV retograde : arc adeb ( e < 1 ) 


Table I. Trajectory Types. 
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most of the trajectory, especially for parabolic and hyperbolic trajectories. This 
cancellation does not occur in the case of patched-conics based on inertial frames as 
described in Appendix A (Wadsworth, 1965). 

The Earth-centered translunar conic trajectory may be cither elliptic or 
hyperbolic. The latter corresponds to shorter (light times and is the only case that needs 
consideration for fast trajectories. For the same reason, the moon-centered conic 
approach trajectory is also modeled as hyperbolic in this study. Both the Earth-centered 
and Moon-centered frames are non-rotating. The motions of the Earth and Moon about 
their common barycenter are approximated by coplanar circular orbits. In general, the 
plane of these circular orbits does not coincide with the three dimensional trajectory of 
the spacecraft. (Wadsworth, 1991) 

B. SOLUTION FOR POINT-MOON INTERCEPT 

The least complicated form of a first order solution is for a point-moon intercept. 
Figure 3 illustrates the relevant parameters of the problem. Appendix B, utilizing 
MathCAD 3.0, numerically solves the first order solution of point-Moon and massless 
spherical Moon intercept. In Figure 3 the following assumptions are made: 

• The Earth-centered cartesian coordinate frame is non-rotating (indicated by 
X,Y,and Z). 

• The Moon-centered cartesian coordinate frame is non-rotating (indicated by x,y,and 
z with the x-axis parallel to the X-axis, and the z-axis tilted at angle i M for the Z- 
axis.) 
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The initial point-Moon intercept solution represents a simple three-dimensional 
model; three-dimensional considerations arc incorporated in the first-order solution solved 
in Appendix B. These considerations will be discussed in section B.2. 

1. Three Dimensional Point-Moon Model 

Figure 3 illustrates the three-dimensional model of the first-order solution at 
point-moon intercept while Figure 4 further illustrates aspects of the three-dimensional 
model. Figure 5 delineates a flowchart of tasks to be accomplished in the three- 
dimensional model of the first-order solution of point-Moon and massless spherical Moon 
intercepts. This flowchart summarizes steps used to construct Appendix B. Symbols and 
abbreviations are defined in the Table of Symbols and Abbreviations. The hyperbolic 
trajectory yielded an energy requirement of 63.2 km 2 /sec 2 while the Apollo mission 
required a minimum energy requirement of 61.8 km 2 /sec\ This results in a 2% increase 
in energy by using the hyperbolic trajectory, and as a consequence, will increase launch 
costs. Considering that the time of flight is reduced by more than 50%, the 2% energy 
increase penalty seems feasible for this mission. 

2. Three Dimensional Considerations 

Figure 6 illustrates the relationship of orbit eccentricity and spacecraft true 
anomaly at moon intercept. Figure 7 is an exploded view of Figure 6 for user 
interpretation which shows eccentricities between 1.0 and 2.0. These diagrams must be 
referenced by the user to insure appropriate inputs for solution in Appendix B. Figure 
8 illustrates a cone sweeping out a circular area of the Earth. The Earth surface within 
the cone defines a zone of inaccessibility for Earth launch sites, given the specified true 
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Figure 3. Three dimensional point moon model. 



Figure 4. Three dimensional model focusing on the point 
Moon. 
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Figure 5. Translunar Trajectory Program Flowchart. 


anomaly at moon intercept. This fact constrains inputs to the program in Appendix B. 
If the launch site is found to be in the restricted cone region a different trajectory must 
be used. 
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Figure 6. Eccentricity versus maximum and minimum true 
anomaly at Moon intercept. 
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Figure 8. Three dimensional considerations. 


C. SOLUTION FOR MASSLESS SPHERICAL MOON INTERCEPT 

The final solution in Appendix B, is the massless spherical Moon intercept. An 
iterative technique is employed to obtain an impact on the lunar surface. This technique 
is illustrated and described in section C.l. Selenographic coordinates are utilized to 
obtain the solution for the massless spherical moon intercept and are described in section 
C.2. 

1. Iterative Technique for Lunar Surface Impact 

Figure 9 illustrates the method to ascertain the lunar surface impact point 
iteratively. The iterative process is defined in the following steps: 
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• Determine point-moon intercept coordinates (R=p EM , f M , <£ M , and t M are inputs). 

• Select trial value of radial distance : R 0 ss p EM -D. 

• Calculate the corresponding flight lime from perigee, t 0 , and separation flight time 

D 0 - 

• Find the next iteration : R, = R 0 + (D 0 - D)cosA 0 . 

• R n *D n -i-R n D n = cosA.,, and R n = R n . 1 + (D n .,-D)cosA. n . 1 (bold print represents vector 
quantities). 

• Continue until n = 10 or | D n .,-D | < 1 km. 


This iterative technique is employed in Appendix B, however, solutions are calculated 


manually by the program user without a terminal loop. 



Figure 9. Iterative technique model. 
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2. Selenographic Coordinates 

Impact points on the near-side of the moon are expressed in selenographic 
coordinates, for which the zero degree longitude and latitude points at the Earth. The 
North, South, East, and West directions indicated on Figure 13 for latitude and longitude, 
are for an observer on the Earth. The lunar librations (Roy, 1988) are neglected in this 
definition. 

a. Euler angles 

Figures 10 and 11 express Euler angle transformations which are 
incorporated in Appendix B to relate selenographic coordinates to Earth-centered 
coordinate systems. The first Euler angle,(0), is zero because this angle is not required 
in this analysis. The second Euler angle rotation,(0), is the angle i M in Appendix B. The 
third Euler angle rotation,(T M ), is the angle <£ M in Appendix B. (Goldstein, 1951) 

b. Lunar Impact 

Figure 12 illustrates a selenographic coordinate solution for lunar impact 
of the spacecraft. Since time reversal is employed in Appendix B, the Moon moves 
backwards in its orbit for the translunar trajectory'. The sign on co M is reversed in 
Appendix B for this reason. 

The selenographic latitude for surface impact is designated 5 
and the selenographic longitude for surface impact is designated a. Figure 13 illustrates 
these latitudes, longitudes and their signs. The solution for a has a factor of it added to 
it to insure that the positive X-axis points to the Earth. Figure 14 displays a locus of 
impact points for varying time offsets (t off ). The axes units are in radians. Varying t off 
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Figure 10. Second Euler angle rotation. 



Figure 11. Third Euler angle rotation. 
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Figure 12. Lunar centered frame. 


SELENOGRAPHIC 

COORDINATES 



Figure 13. Selenographic coordinates. 
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is equivalent to varying the input value of the lunar sweep angle (<£ M ). By providing an 
offset from the point-Moon intercept declination (<5 M Q(1 ), the user can target from North 
to South (see Figure 13). This factor can be combined with (t off ) to produce an arbitrary 
location of the impact point on the lunar near side. 
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Figure 14. Locus of lunar impact points. 


D. PATCHED CONIC APPROXIMATION 

The method of calculation for the massless spherical moon intercept model can also 
be used for sphere-of-influcnce calculations if the radius of the sphere-of-influcnce is 
substituted for the radius of the Moon. The patched conic lunar impact point, using a 
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finite lunar mass model, can then be approximated by an iteration process identical to the 
process used in Appendix B (Bate et al., 1971). 

From the final patched conic iteration, the impact point for the Moon-centered 
conic trajectory can be calculated. This calculation is made using the same method as 
Appendix B. 
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MI. APPLICATIONS 


Fast Earth-return trajectories have numerous applications, such as for requirements 
in evacuating a lunar base in an emergency or the establishment of a lunar-based 
strategic deterrent force. This chapter will discuss such applications, their development, 
and technical feasibility. 

A. STRATEGIC PEACE INITIATIVE CONCEPT 

The ABM treaty prohibits nuclear weapons in space. If we look at a future world, 
this treaty might possibly be altered or not exist. The Strategic Peace Initiative is a 
hypothetical plan to place nuclear weapons on the Moon and to eliminate them from the 
Earth. The reasoning for this replacement is to prolong the strike capability and reaction 
time that now exists, thereby increasing stability. 

1. Background Information 

This section describes a hypothetical Strategic Peace Initiative, a plan to 
eliminate the risk of global nuclear war and simultaneously promote peaceful uses of 
space. The SPI concept could be a step toward reducing present nuclear confrontation 
as well as "denuclearizing" the Earth. (Wadsworth, 1988) 

Rather than using our current posture of SDI, the SPI concept can encourage 
peaceful exploitation of space (Gray, 1985). SPI resources could provide the logistics 
support for a new era of nonmilitary space activities, such as manned space exploration, 
servicing geostationary communication satellites, safe disposal of radioactive waste, clean- 
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up of hazardous space-debris, and eventually, deflection of Earth-threatening asteroids 
or comets. 

Although the basing of deterrent strike forces in space, particularly on the 
Moon, has already been proposed; this proposal is quite unique. The novelty of an SPI 
concept is that by a new treaty, the United States and Soviet Union would share the 
"high ground" of space. As with both the INF and Start treaties, success of the SPI treaty 
would depend on whether it serves the mutual interests of parties involved. (Wadsworth, 
1988) 

The free world’s offensive deterrent policies have prevented a global nuclear 
confrontation to date. Unfortunately, this policy is dangerously unstable because of the 
critically short decision time for launching a counterstrike to a perceived nuclear attack. 
SPI would eliminate or greatly reduce the danger of accidental nuclear war by increasing 
decision time from minutes, in the case of an SLBM attack, to a more comfortable two 
days. 

The primary justification for SPI is survival in a nuclear age, consequently the 
emphasis in this proposal is on deterrence. A secondary, economically-compelling 
justification is that much of the SPI investment in global security could also serve the 
civilian space program. They both require similar logistics elements: 

• A heavy-lift launch vehicle 

• A space tug (space/lunar transfer vehicle) 

• Earth/Moon space communications and tracking networks 

• A manned lunar base 
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2 . 


SPI Architecture 


The minimum-energy transearlh strike trajectory has a 4.7 day flight time. 
Doubling the injection energy yields a three day trajectory. High-energy strike trajectories, 
taking two days or less, still permit nuclear strike weapons to escape the Earth-Moon 
system if aborted. Appendix B provides such a time of flight (approximately 1.9 days). 
Shortening the flight lime from 4.7 days to 2 days results in only a two percent increase 
in transearth trajectory energy, which is equivalent to providing the additional velocity 
increment typical of an 1CBM. Atmospheric re-entry velocity does not increase 
significantly between 4.7 day and 2 day trajectories, which simplifies re-entry vehicle 
design. To minimize vulnerability to intercept, the trajectories must be of the direct- 
ascent type without lunar parking orbits, in contrast to the more economical supply 
shuttle trajectories which must also be utilized. (Wadsworth, 1988) 

The SPI communications networks and surveillance systems would preferably 
be located at the natural synodic satellite points, designated LI through L5, which are 
stationary with respect to the Earth-Moon system and would minimize station-keeping 
fuel. These points are therefore an ideal location for placing these communications 
satellite systems. (Farquar. 1970) 

The two opposing lunar base zones, each approximately one hundred 
kilometers in diameter would contain the strike force missile silos and manned 
command/support complexes. These zones would be located on approximately opposite 
sides of the moon. They might occupy polar sites to take advantage of full-time sunlight 


for solar power. Determining optimum zone locations requires analysis if the optimum 
strike trajectory missions for both bases. (Wadsworth, 1988) 

Similar to the deployment of the stages of SDI, SPI would be implemented 
in stages and be subject to treaty compliance verification (possibly a United Nations role) 
at each stage. Both superpowers would deploy the following five stages: 


• NON-NUCLEAR DEPLOYMENT : Lunar bases, communications, and 
surveillance systems deployed. 

• FLIGHT TESTS : Test strike force with dummy warheads to demonstrate accuracy 
and reliability. 

• NUCLEAR DEPLOYMENT : Equip lunar bases with nuclear missiles and phase 
out Earth-based strategic nuclear forces. 

• VERIFICATION : Surveillance of lunar bases and Earth-Moon space; inspection 
and monitoring of Earth-based nuclear device manufacturing and test facilities; and 
inspection of translunar cargos. 

• STRATEGIC REDUCTION EVOLUTION : Bilateral reduction of strategic 
nuclear forces to minimum required for asteroid deflections; increase in cooperative 
nonmilitary space ventures by utilization of lunar bases. 


B. SPI TECHNICAL FEASIBILITY 

Initial deployment of SPI relies on technology that has or will soon be proven by 
other programs. In fact, SPI can virtually be used as a spin-off of developed SDI 
technology. Programs such as TDRSS, NASA space station, OMV, NASP, HEDI, and 
the APOLLO lunar programs are examples of some of these programs. 

Challenges in technology' do exist. These challenges include a precision guided re¬ 
entry' vehicle, developing "smart" surveillance techniques, and prevention of adverse 


physiological effects of long-term residency in low gravity. At an interval of deployment 
often to twenty years, the fully deployed SP1 deterrent force system probably would cost 
comparably to the proposed SD1 system. (Wadsworth, 1988) 

Under the SP1 treaty, a sovereign lunar base zone would be assigned to each 
superpower for basing its nuclear strike force. These zones would be located on 
approximately opposite sides of the moon. The optimum size and location of the zones 
would be a compromise among the following six basic requirements: 


• Minimize Surveillance Cost: Requires minimizing zone size; locale bases on the 
Moon’s visible side to take advantage of surveillance redundancy provided by Earth- 
based sensors to back up the spaced-based sensors. 

• Maximize Nuclear Survivability: Requires dispersal of hardened missile silos within 
a circular area of at least one hundred kilometers diameter. 

• Minimize Base Vulnerability: Requires maximizing separation between bases to 
deter pre-emptive strikes and reduce vulnerability of transearth strike forces to 
interception. 

• Minimize Transportation Costs: Requires selection of base zone locations to 
minimize fuel cost for both supply shuttles which utilize parking orbits and strike 
forces which require direct-ascent trajectories. (This issue will be difficult to resolve, 
because if both strike forces were to take advantage of the Moon’s orbital velocity, 
then locations near the lunar poles will be favored. This induces a strategic 
problem, because the energy requirements for striking Northern-hemisphere Earth 
targets will differ between the polar bases.). 

• Avoid Base Overflight: Requires designing supply shuttle parking orbits not to 
overfly opposing base. 

• Minimize Communication Costs : Requires designing base communications on the 
near-side of the Moon. Far-side Moon bases require more communication system 
elements because its orientation with the Earth docs not meet full-time line-of-sight 
communications link requirements. 


C. SPACE SUPPORT ACTIVITIES 


A large variety of space support programs can be served by SPI. For example, 
lunar-based astronomy would be free of restriction imposed by Earth’s atmosphere. A 
lunar base could be used as a staging platform for a manned expedition to Mars, taking 
advantage of a smaller round trip propulsion requirement for liftoff from the Moon, 
rather than from the Earth. The servicing of geostationary satellites could be 
accomplished with SPI resources. Transport of concentrated nuclear waste to a safe 
dumping site in the Moon would be unhampered by groundwater problems as on the 
Earth (Rosen, 1981). SPI resources could be used for debris removal from low-Earth 
orbits. This is an increasingly serious hazard to space missions. Also, safe disposal of 
space junk or dead satellites in higher orbits could be an SPI priority. 

Using the SPI surveillance network, the nuclear arsenal of the Moon can be used 
to deflect Earth threatening asteroids or comets. Collision with a ten-kilometer diameter 
asteroid has the explosive force of ten million one-megaton nuclear bombs. With an 
explosion of such magnitude, the impact dust cloud will plunge the world into a deep 
freeze similar to "nuclear winter” but much more devastating. Fossil records indicate that 
such an asteroid could account for extinction of the dinosaurs and cause for the Ice Age 
(Malove, 1985). 

D. SUMMARY 

The Strategic Peace Initiative Concept is a revolutionary plan for strategic nuclear 
deterrence. It is a plan which provides double value; a solution to accidental global 
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nuclear war by increasing reaction time to at least two clays, and promotion of peaceful 
uses of space by developing permanently manned lunar bases and their infrastructures. 

The uniqueness of SP1 is that it requires bilateral treaty agreement between the 
United States and the Soviet Union. Although the politics of such a treaty requires 
extensive thought, verification of the treaty can be accomplished. By looking to the 
future of peace, perhaps we can be at peace with ourselves and with future generations. 
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IV. CONCLUSIONS AND RECOMMENDATIONS 


This final chapter summarizes the findings if this thesis. It also recommends 
research for further work concerning this topic. 

A. FINDINGS 

This thesis utilizes the selection of an efficient and operationally meaningful set of 
inputs and terminal parameters, to lead to a transearth or translunar trajectory targeting 
solution, which rapidly converges in an iteration process. 

B. RECOMMENDATIONS 

This thesis lays the groundwork for an accurate patched conic method of designing 
fast transearth trajectories. Recommendations include the finalization of the patched 
conic approximation and the iterated impact point for the Moon-centered conic 
trajectory. A development of the SPI concept to include a cost/mission analysis would 
also be very useful. 
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APPENDIX A. PA TCI I ED-CON IC BASED ON INERTIAL FRAMES 


The patched conic trajectory can be formulated in terms of Earth-centered and 
Moon-centered coordinate frames which tire inertial, being attached to a Fictitious Earth 
and Moon which travel at constant velocity with respect to the barycentric frame. The 
Earth-centered conic trajectory approximation implies the lunar gravitational force is 
ignored so that the Earth is not accelerating with respect to the barycentric frame. 
Consequently, after the instant the spacecraft is launched, the center of gravitational 
force, the Earth, is replaced by a fictitious Earth which is imagined as moving in a straight 
line at constant velocity in the barycentric frame. In the meantime, the massless Moon 
continues on its circular orbit in the barycentric frame. In terms of the translating, Earth- 
centered frame, the lunar orbit appears as a compound curve, rather than a circle. This 
fact must be taken into account in targeting the conic trajectory for lunar intercept or 
crossing of the lunar-sphere-of-influence. The same considerations apply to the Moon- 
centered conic phase where the Moon is approximated by a fictitious Moon moving in a 
straight line at constant velocity in the barycentric frame. For greatest accuracy, the line 
is taken tangent to the Moon's circular orbit in the barycentric frame at the time of 
intercept. This requires an iterative solution. (Wadsworth, 1991) 
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APPENDIX 15. LUNAR TRAJECTORY PROGRAM 


input nmmm 


S/l' declination at caitli !aunc!i point 

0, n • 28.3 

180 

(input in degrees) 

Oy = 0.493928178314395 

radiaiiS 

true nncmely of s/c el moor intercept 

*164 

M 180 

(input in degrees) 

f M = 2.8623399732707 

radians 


inclination o! the moon orlil plane 
(input in degrees: 

i M • — * 20 

M 180 

20.5 > i > 10.5) 

M 

i M = 0.349065850398866 

radians 

eccentricity of the transkna? trajectory 

(eccen > = ecccn ) 
min 

ccccn : 1.01 


geocefilric altitude of launch point 
(altitude in kilometers) 


He - i 


lunar weep anpje el tnocri intercept 



(input iti degrees) 


<f, M = 2.96705972839036 
radians 


posigrede or retrograde orbit determination 
posdgrade = +1 

retrograde = -1 
“E obt = 1 


(not iterated in program ) 


5 [4 oft 0 


ramie OUTPUT: 

0 M ■= asin[sin[i M ]‘sin[* M ]] 

6 m = 0.059426145347819 
(spacecraft declination cl moon intercept) 

_ cos Um] 

“ M1 " [cos[o M j 
a mi =-0.986549223316463 



«M2 = 0.163464460888851 


: = 2*atan 


[ 1- “M1 ] 

"M2 


u M = 2.97739131857372 
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INPUT CONSTANTS: 


pem = 384400 km : earth-moot meat distance 
r£ := 6378 km : earth mean rad its 

r M := 1738 km : moot meat radios 

cj E := 7.292*1 o' 5 rad/seo: earth atgtlar relation rale 

COM := 1.525-10~ 4 *— 

M 180 

<o M = 2.661627109291352-10 -6 ni/sCC : 
moot mean orbital rolalional rale 

2 3 

k£ := 398600 km /sec : earth gravitational constant 
2 3 

k M : = 4903 km /see : moot gravilaliotal constant 

OUTPUT VARIABLES and DRMUUS: 

<o M := -i*w M (dte to the moot moving backwards it its 

orbit lor Uc translator trajectory) 

R E1 ;= r E +h E 

R ei = 6379-10 3 km : (s/c radial distance at earth launch) 
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R M1 PEM 


r m1 = 3.844M0 5 km : (s/c radial distance at mooi intercept) 


r 


R 


Ml 


Rf 


El 


(radial distance ratio} 


r = 60.26022887599937 


fir = acos 


r"£ 1 +cccerTcos 


I'm]] 


" 1 


eccen 


f E = 0.726288833059664 


f EM - f M~ f E 


f EM = 2.136051140211037 

si £ n fEM ;= sin [ f EM] 
sign fEM = 0.84445214187218 

Sfem ■“ 1 


NOTE: for S {cm 

if sign is positive then use 1-1 
if sign is negative use -1 


“EMI 


[ cos [ ^ EM ] ~ sin [ ^ E ] * s ‘ n [ ^ M ]) 
[cos[5 e ]-cos[6 m ]] 


“EMI =-0-641452764362017 
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“ EM2 


F 


EMI 


* s fcm* a E OBT 


a EM2 = 0.767162532383019 


u EM 2*alan 


[ 1-a EMl ] 

“EM2 


NOTE, 

if a EM1 = -1 then 

is *'aE_OBT 

if “ EMI = ^ 11)611 

is 0 


a EM = 2367186781944015 


“ E “ M ' “ EM 


a E = 0.710204536629705 (s/c rigtl axeosics al earth laonch) 


$El : = -ir + fM-«M 
6 E1 = -0338678825666912 

*E2 := *~ f M -5 M 
5 E2 = 0.219826534971273 


f M_MAXl = sin [ 5 E]‘ sin [ 5 M] + cos [“EM]* cos [ i E]* cos [ 5 M] 


f M_MAXl =-0-53563101113312 


f M_MAX2 - J 1_I M_MAX1 2 ] 


( M_MAX2 = 0.84445214187218 


“EM 


“EM 
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f M_MAX 


2*atan 


[ 1-{ M_MAX11 

f M_MAX2 


f M_MAX 

= 2.136051140211037 


If: 

5 El <5 E <5 E2 


then 

f M <f M_MAX 

is constrained 


oUenn9e cse 


A a := asin 


tan[d E ]*tan[5 M ]- 


COS 

f M 

l 

M 

[ 5 e 

*cos| 

:*««] 

i] 


Aa = 1.570796326794897 - 0.4963212320603151 


Coofllraicl: 

aE MIN := a M~y Aa 

aE MIN =-0-164201335016073 + 0.496321232060315i 


a E m«sl be greater tfean 


aE MlN 


eccen 


mm ' 



r 




_ . 


_ 


1“ 

2 

' cos [ 1 EM ] + 

1 


r 

2 



1 


r 


ccccn = 0.974647175723671 


a E 

as calculated 
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energy m i n 



1 + ecccn 

min ]] 

[ 

;2-r M i ] 

i 


( = minifnvm specific energy 
at buck} 


energy min = -1.023795999276086 


sin 

l*«] 

M 

5 E 

J*cos| 

[ { em] 


cos 

1 § E. 

l‘ S,n l 

. f EMj 

1 


j9 E1 = 0.421410460300705 


0E2 : " 


COS 

[•■*: 

]‘ sin l 

[“em] 

sinj 

. { em] 

1 


P E2 = 0.906870014913464 


/?£ 2"atari 


1 ~^E1 
0E2 


flight path azirioUi at lenicti: 

TC>P>~% 


= 1.135796261808675 


\\ cosj^6 £ j*sin^£ j 


ij = 0.798478510707822 


>2 J 1_ 'l 2 

i 2 = 0.602023311789352 
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i 2 # atan 


(InansloDar trajectory iidinaUoi: 
0 < = i < = ir) 


l-il 


i = 0.646032654504068 


aE fll 


A El 
'2 


oEqj = 0.699990269559789 


aE Q2 


tan[a E ] 

tan(i) 


oEq2 = 0.714152380463452 


ocEq 


2*atan 


aE Q2 


uEq = 0.795412455434794 

O := a E -aE fi (traidvDar trajectory loogitvde- 

of-aaeDdinf-Dode) 


£i = -0.085207918805089 


CEl cos [ ^ E ]*cos[aE q ] 


C E1 = 0.616325580024169 
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c E2 := 


sin[ 5 E ] 


C £2 = 0.787491447196649 


C i 7 •- 2“atari 


1-C 


El 


' E2 


Cr = 0.906728172342672 


(s*eep angle at laoccfc: 
Sr>c > = 0) 

e 


Rp 


R E 1 1 +eccen*cos£ f E ]] 
1 + eccen 


(perigee distance) 


Rp = 537010652841956-10 


V P = 


jk.£ ■( 1 +cccen) 


R, 


(perigee velocity) 


Vp = 11.99319583652338 


| 2>k E 


( 1 - ccccn) *k.£ 

J1 r ei j 


r p I 


(s/c velocity el laoidi) 


V E = 1131107362995628 
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Ra : = 


R p*( 1 -eccen) 
1 - ccccn 


(apogee distance; valid only lor eccen ( 1 
ignore for eccen > = l) 


R a = 5.57010652841956 *10 3 


^asy ;r ^~atan 




eccen - 1 


fasy = 3.000756780023376 


(hypertwlic esytnplolic tree anomaly 
ball-angle; valid only lor eccen3 1 
ignore for eccen < = 1) 


D : =2*f aS y-x (deflection angle) 


D = 2.859920906456958 


d asy ’ — ^ p 


1 + 


eccen 


= 1.108506348725081 *10 4 


(directrix-to-locos didance or axial olsel 
of asymptotic from locos; valid only lor 
eccen 3 1 ignore lor eccen < = 1) 


7E1 := 




Y£1 = 0.934111286367402 


VE2 


eccen*y E1 * sin^ f £ j 
1 + ccccn *cos[ f £ 1 


y E2 = 0356981938872315 
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y£ := 2‘atan 


yE2 


(digit path angle at leincti; 


->tf >— 
2 h 2 


y E = 0365034946775507 


M : 


K £ 

( 1 - ccccn) 

| R Ml 

R P J 


(s/c velocity at moon intercept] 


V M = 1.670175762283258 


y mi := 


r p‘ v p 

r mi-v m 


y M l =0.104052584194718 


f Ml = C03 [ { M] 


[ M2 sin [ f M] 


7M2 : ~ 


cccen-y M l 
l+[ ccccn'f mi ] 


y m =0.99457179716821 
* - 7M1 


y M := 2-atan 


y M2 


(digit path angle at moon iilercept) 


y M = 1.466555060009538 
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Pm := ‘2*[ C0S [ a M' Q ]] 


0 M1 =-0.600145987662734 


0M2 


'1 



&M2 = 0.799890488437212 


*— 2'atan 


1 ~^M1 
&M2 


(fligkl path azimntk at moon intercept) 


0 M = 2.214479932657055 


X M := R M1 ‘ C0S [ d M]‘ C03 [ a M] 

(RECTANCUUR COORDINATES OP 


THE KOON AT POINT-MOON 

X M = "3.785601002578928 MO 5 

INTERCEPT) 

y M - R Ml ,cos [ 5 M]* 3in [“M] 


Y m = 6.272482025241596-10 4 


Z M := R M1 * sin [ 5 M] 


Z M =2.282996752157727 *10 4 



u Ce~ 1 e 


(traistonar Irajedory aigcmeil- 


w = 0.180439339283008 of- perigee} 
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c M f M +w 


Cj^ = 3.042779312553709 (sweep angle al moon intercept) 


k M' 


M 


COS 


*[[7M2“ sin [ a M]‘ i 2‘ C0S [ c M“7M]]‘ C0S [ a M" Q 


[ 5 m] 


X M - = -1.665737958076395 


M’ 


VM 

cos[ 6 m ] 


*[[7M2~ sin [ 5 M]‘»2‘ cos [ c M"7M]]* sin [ a M" Q 

* C os[a M -Q]] 


Y M - = 0.048513058379955 
Z M’ := v M*'2‘ cos [ C M"7M] 


Z M - = -0.00545766973131 


r M’ 


r mi 


(redial Telocity al poinl-moon intercept) 

R M ’ = 1.648023680688986 
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2-ir‘R 

t •=_ L 

P . 5 


1.5 


(circclor orbil period for radios R) 

P 


T p = 4.137199589247244 * 10 


1 - cos 


‘El 


t'E] 

in [ f E ] 


f E1 = 037999703209726 


f Ml := 


1 "COS 


Cm] 


sin 


Cm] 


f M1 = 7.115369722384208 


E h i -= In 


1 + 

> 

ccccn - 1 , 

1 FI 

1 ♦ ccccn 

'■J 

ccccn - 1 _ 

El 

1 + ccccn 


E hl = 0.05361869041727 


E h2 := ln 


. (ccccn “ 1 . 

1 + h-*f Ml 

1 + ccccn 


ccccn - 1 
1 + ccccn 


•f 


Ml 


E h 2 = 1.103630544045282 




















-shl r 


ecccn ~1 . 

1“ -*f f i 

1 + eccen 


-1 

/ 

- 

[ ,4 J 


ccccn -1 
1 +cccen 


ri 


*f£l 


E shl =0.053644386077804 


^sh2 


1- 


Icccen - 1 
^ 1 4- ccccn 



1 + 


ecccn - 1 
1 + cccen 


-1 


Ml 


E sh2 = 134171391602585 


_[ccccn-E shl ]-E hl 

* El l5 

2‘ir*(ccccn - 1) 

t E1 =0.089467283524084 


t 



ccccn ‘ E 5h2]~ E h2 
2*ir*(eccen- 1) L5 


t M l = 40.02754953820076 


1 E - Tp^El 


(earth laoucl Right lime relative 
to perigee passage) 


t E = 370.144008646907 
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1 M := T p M Ml 


(poiil moan intercept fbgfat tine 
relative la perigee passage) 


t M = 1.656019615080179 ■ 10 S 

* EM •= t m _ t e (peiil-moon Iratishnar trajectory 

flight lime) 

*EM = 1.65231817499371 *10 5 


INPUTS: 

R 0 ; = 382986.2976625 (aiy radial distance) 


OUTPUTS: 

v jt 2 ‘ k E ]_[< 1 _ccccn )‘ k E] 

°"J Ro " R P 


V 0 = 1.6724659380983 


R. 


_ [R p *(l ♦«*«!)] 


Rr 


R x = 0.029233197611653 


fo := 2'atan 


ccccn + 1 - R x j 

l5' 

ccccn - 1 + R x j 



f 0 = 2.861954051456415 
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c0 - w + fo 

C o = 3.042393390739423 


5q asin^sin(i)*sin^Co]] 


5 0 = 0.059657751937721 


- H-B 

[i7n(ij-Un[c M ]] 


Aa M1 =-0.996881642139325 


Aa ui .= —i- i 

Mli tan(i) 


Aa m2 = 0.078911289214964 


Aa m := 2’atari 


[* ~^**Ml ] 

Aa M2 


Aa M = 3.062599237378809 


1*4*01 

A/y 1 :=_ L J 

[sin(i)*tan[C 0 ]] 


Aa t = -0.996857191524672 
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A a 2 


lan[5 0 | 
tan(i) 


Act 2 = 0.079219566431173 


Act 0 


2*alan 


[1-AajJ' 

t±a 2 


Aa 0 = 3.062289992046826 


“j := “M _ Aa M 

a l = -0.085207918805089 


a 0 := Aa 0 + a t 


a 0 = 2.977082073241737 


[ R P* V P 

y, "pVV^ 


Yi = 0.104293660037493 


i 

i 


ccccn'y j *sin 

[*o] 

I] 

[ 1 +cccen*cos| 

N 

!] 


y 2 = 0.994546546158593 
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Yq := 2*atan 


i'-ni 

>2 


y 0 = 1.466312665340157 
P\ - i2*cos[a 0 -Q] 

= -0.600131267822715 



p 2 ~ 0*799901532303508 


Po : = 


2*atan 


f‘-M 


p 0 = 2.214461530464985 


Xq '■= Ro*cos[6 0 ]* cos [ a o] 

X o = -3.771433240671579*10* 


(RECTANCVUR COORDINATES OP 
THE SPACECRAFT ON ITS TRAJECTORY) 


Yo : = Ro ,cos [ 5 o]* sin [ a o] 


Y o = 6.260990806241751 MO 4 

Zq : = Ro*sin[6o] 

Z 0 = 2.283455104046146 *10 4 
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1 - cos[ f o ] ] 

sinpo] 


f mo = 7.105421102050813 


h3 


:= In 


1 + 

((ccccn - 1) 

* f M0 


^J( 1 + ccccn) 

1- 

(ccccn ' 1) 

* f M0 

1 + ccccn) 


E h3 = 1.101755458701495 



r 

--- - 


"I 

r r 

^sh3 * 

1~ 

)(ccccn - 1) 
/\J( 1 + ccccn) 

‘ f M0 

- 

1 + 

. 



. 


J (eccen - 1) 
(1 ♦ ccccn) 


f M0 


E sh3 = 1-338578545032142 


*M0 


[[ccccn»E sh3 ]-E h3 ] 
[2*x*(cccen~ l) 1-5 ] 


t M0 = 39.8219787493874 


*0 “ V‘MO 

t 0 = 1.64751474124978 *10 5 
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'off - 0 


(lime ofsel awl) 


Atq := 'M _t 0 +t of! 


Ar 0 = 850.4873830398719 


♦ Ml := C0S [ a M ]‘ C0S [^M ] 


♦ Ml ="0-984807753012208 


♦ M2 ; " 


sin [ ^ M ] 

sin[i M ] 


♦ M2 = 0.17364817766693 


♦m := 2,aUn 


[ '"♦Ml ] 

♦ M2 


= 2.96705972839036 


♦ MO ;= ♦M _w M‘ Ar 0 


* M0 = 2.%9323408665169 
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X M0 - PEM‘ cos [*MO ] 


(RECTANGULAR COORDINATES 0»’ THE 


X M0 =-3.787102316833942 MO 5 


YmO = PEM * cos [ 'M ]* s ' n [ ^MO ] 


Y M0 = 6.191940094679928 *10 4 


Z M0 

Z M0 = 2.253681886822944-10 4 


NOON AT TIUC ( } 
0 
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D M 0 ; = 4[X 0 -X M of + [Y 0 -Y M o] 2+ [Zo-ZMof] 


DISTANCE BETWEEN THE SPACECRAFT AND THE LUNAR CENTER AT TIRE l 

0 

D M0 = 1.738000000264939 *10 3 km 

' THE COAL IS 10 DRIVE THIS DISTANCE TO EQUAL THE RADIUS 
OF THE MOON FOR AN IMPACT ON TIE SURFACE OF THE MOON 


: = D M0~ r M 

error = 2.649394446052611*10 


[[XmO'Xq] 

|- x o*l 

[Ymo-Yq]* 

y o + [ z mo-Zq] 

N 

O 

[ r o‘ d mo] 

1 


Apo = 0.81263979792074 ^ cos [ x n] ^ 


^next “ *DMO‘ApO 


R nc „ = 3.843986656315015 MO* <«*»*« l»*«l »«*« •' 

R 0 ) 
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X := X M -X M0 


X = 150.1314255014295 


Y =Y M -Y M0 


Y = 805.4193056166769 


Z z M _z M0 


Z = 293.148653347831 


INPUTS: 

♦ := o 


( EULER ANGIE INPUTS) 


e := i M 


= *mo 


♦ M = 2.969323408665169 
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OUTPUTS: 


(THE MATRIX STRUCTURE IS NUMBERED AS FOLLOWS:) 


A i B i c i 
A 2 b 2 C 2 
a 3 b 3 c 3 


( ROTATION AND TRANSLATION ) 


A, := cos[> M ]‘cos(^)-cos(0)*sin(*)*sin[¥ M ] 
A { = -0.985198313432347 

B £ := -sinj^ ♦ J*cos(^) -cos(0)*sin(^)*cos[^ ♦ m ] 

B! = -0.171418444777856 

C j := sin(0) 'sin(^) 

C, = 0 

A 2 := cos[ + M ]*sin(^) + cos(0)*co5(^)*sin[^ M j 
A 2 = 0.161080647624348 

B 2 := -sin[i' M ]‘sin(^)+cos(0)*cos(^)*cos[>P M j 

B 2 = -0.925783585143099 
C 2 := -sin(0)*cos($) 

C 2 = -0342020143325669 
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A 3 •= sin(0)'sin[* M ] 


A 3 = 0.058628561051585 


B 3 := sin(0) *cos[ ♦ M ] 


B 3 = -0336957668364338 


C 3 := cos($) 


C 3 = 0.93%92620785908 


A := 


‘1 


B 1 C, 
a 2 b 2 c 2 
A3 B 3 C 3 


B := 


X 

Y 

Z 


C := A*B 


C = 


-285.9729519600736 

-821.7234494824876 

12.87940451383525 
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C 0 = -285.9729519600736 


Cj = -821.7234494824876 


C 2 = 12.87940451383525 


^ lun C 0 

Y| un : = C, 

^lun = C? 


6 := atan 


'lun 


^^lun + ^lun 


6 = 0.014801756974441 


SEIENOCRAPHJC LATITUDE FOR SURFACE 
IMPACT ON THE NOON 


a 



+ * 


a = 1.235890035212915 


SELEN0CRAPH1C LONGITUDE FOR SURFACE 
IMPACT ON THE NOON 
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